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The metaphase spindle is a dynamic bipolar structure 
crucial for proper chromosome segregation, but 
how microtubules (MTs) are organized within the 
bipolar architecture remains controversial. To explore MT 
organization along the pole-to-pole axis, we simulated 
meiotic spindle assembly in two dimensions using dy- 
namic MTs, a MT cross-linking force, and a kinesin-5— like 
motor. The bipolar structures that form consist of anti- 
parallel fluxing MTs, but spindle pole formation requires 
the addition of a NuMA-like minus-end cross-linker and 



directed transport of MT depolymerization activity to- 
ward minus ends. Dynamic instability and minus-end 
depolymerization generate realistic MT lifetimes and a 
truncated exponential MT length distribution. Keeping 
the number of MTs in the simulation constant, we ex- 
plored the influence of two different MT nucleation path- 
ways on spindle organization. When nucleation occurs 
throughout the spindle, the simulation quantitatively re- 
produces features of meiotic spindles assembled in Xeno- 
pus egg extracts. 



Introduction 

The spindle is composed of microtubules and hundreds of 
other factors, but how these components interact to generate a 
dynamic, steady-state structure is unclear. Spindles assemble in 
Xenopus egg extracts through a self-organization pathway in 
which a gradient of RanGTP surrounding chromatin promotes 
microtubule (MT) nucleation and organization of MTs by 
motor proteins into an antiparallel, bipolar array. Although 
the spindle appears static, it is actually a steady-state structure 
composed of short-lived MTs. MTs align predominantly along 
the pole-to-pole axis and slide poleward through the action of 
kinesin-5, a plus end-directed tetrameric motor (Walczak and 
Heald, 2008). At the spindle pole, localized MT destabilization 
results in depolymerization from MT minus-ends. During the 
metaphase steady state, sliding and pole-localized disassembly 
are coordinated to form a spindle of constant length exhibiting 
"flux" (Dumont and Mitchison, 2009). Yet, how such a balance 
between MT assembly, sliding, and disassembly is achieved is 
not understood. 

The contributions of various mechanisms to spindle archi- 
tecture remain controversial, in part because the high density 
of polymers in the spindle prevents current optical imaging 
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techniques from resolving individual MTs. Growing MT plus- 
ends have been observed throughout the spindle, but indirect 
measurements of plus-end dynamics within the spindle disagree 
as to whether MTs near the chromatin are stabilized relative to 
distal MTs (Wilde et al., 2001; Athale et al., 2008; Needleman 
et al., 2010). Moreover, the spatial distribution of MT minus- 
ends remains elusive, and speckle microscopy experiments have 
led to contradictory interpretations of the MT length distribu- 
tion ranging from a truncated normal distribution with a mean 
MT length of ^20 urn (G. Yang et al., 2007) to an exponential 
distribution with a mean of ^3 um (Verde et al., 1992; Needleman 
et al., 2010). Nucleation events cannot yet be observed, and 
multiple nucleation pathways have been proposed to contribute 
to spindle self-assembly. A gradient of RanGTP surrounding 
chromatin may induce nucleation near the midzone (Kalab 
et al., 2002; Caudron et al, 2005), while nucleation along existing 
MTs (amplification) may rapidly increase MT density within 
the spindle (Janson et al., 2005; Mahoney et al., 2006; Clausen 
and Ribbeck, 2007; Goshima et al., 2008; Zhu et al., 2008). 

Computer simulations are well suited to examine the con- 
tributions and cooperativity of these mechanisms in spindle 
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formation because a small number of well-characterized com- 
ponents appear responsible for the bulk of MT self-organization. 
We aimed to simulate important general features of Xenopus 
meiotic spindles, including pole formation, MT dynamics and nu- 
cleation, and spindle length regulation, and then to use this model 
to elucidate MT organization along the pole-to-pole axis. 

Previous models of spindle assembly have not fully incor- 
porated MT plus-end dynamics. In the most recent model of a 
Xenopus spindle, Burbank et al.'s slide-and-cluster model, MTs 
turned over with short lifetimes, but all MTs were of constant 
length and plus ends were nondynamic (Burbank et al., 2007). 
Plus-end growth is difficult to include in a model aiming to pro- 
duce bipolar structures because growth is significantly faster 
than MT sliding. The disparity results in each MT plus-end extend- 
ing three times farther than the minus end has slid poleward. 
If all MTs were nucleated at the midzone, as in the slide-and- 
cluster model, MT structures would be mainly astral: plus ends 
would grow far beyond the region of antiparallel overlap. In this 
work, we explore how dispersed nucleation of dynamic MTs 
resolves this issue to generate bipolar MT organization. 

The mechanism aligning MT minus-ends to form a spin- 
dle pole also remains unclear. The slide-and-cluster model 
aligned minus ends using a dynein-like "cross-linking" motor 
similar to the oligomeric motors shown to cluster minus ends in 
asters (Nedelec et al., 1997). Yet, other work suggests that a 
simple balance of forces between kinesin-5 and kinesin-14/ 
dynein motors does not readily generate MT structures of stable 
length (Badoual et al., 2002; Hentrich and Surrey, 2010). The 
necessary contribution of MT-destabilizing enzymes to pole 
formation (Gadde and Heald, 2004) has not yet been explored 
theoretically. We show here that minus end-directed transport 
of depolymerization activity can produce spindle poles without 
the help of clustering motors. 

Results and discussion 

Simulation of a meta phase meiotic spindle 

We began by populating a simulation with 200 MTs to create a 
two-dimensional (2D) spindle slice of ^5 um in width, on the 
order of internal bundles observed in Xenopus spindles (Mitchison 
et al., 2004). Limiting the model to 2D made the simulation 
computationally tractable while still allowing us to address lon- 
gitudinal MT organization in the spindle, which was robust 
against ±50% variation in MT number (Fig. SI E). To further 
reduce the computation time, MTs were initialized in an anti- 
parallel array (Fig. 1 A, color indicates MT orientation), but 
other initialization conditions produced similar steady-state 
structures. MT minus-ends were inert (Fig. 1 A, white segments), 
whereas plus ends underwent dynamic instability, alternating 
between growth and shrinkage states with parameters in the 
"bounded" regime based on measurements in Xenopus egg ex- 
tracts (Table I; Dogterom and Leibler, 1993; Wilde et al., 2001; 
Brown et al., 2007). Once MTs depolymerized to a minimum 
length (40 nm), they were removed from the simulation. To mimic 
chromatin-driven spindle assembly, we delineated a zone of 
RanGTP activity extending 20 um in length and 5 um in width 
in which MT plus-ends were more stable, incorporating the 



measured downstream effects of RanGTP on MT catastrophe 
and rescue rates (Fig. 1 A, yellow rectangle; Wilde et al., 2001; 
Athale et al., 2008). The number of MTs was kept constant 
through two different mechanisms of nucleation. The first mim- 
icked chromatin-mediated MT nucleation, which is known 
to occur in high concentrations of RanGTP (Kalab et al., 2002; 
Caudron et al., 2005). In the second mechanism, new MTs were 
nucleated along the side of randomly chosen existing micro- 
tubules, modeled after the MT amplification reported to occur 
within spindles (Mahoney et al., 2006; Goshima et al., 2008; 
Zhu et al., 2008). 

In the spindle, many proteins cross-link adjacent MTs, re- 
sulting in an average MT spacing of 50 nm (Mitchison et al., 
2004). Rather than explicitly modeling individual cross-linking 
proteins, we specified that MTs closer than 50 nm repelled each 
other, whereas MTs between 50 and 90 nm apart attracted each 
other (Fig. 1 B). Under these conditions, MTs exhibited Brown- 
ian motion and aligned without crossing one another, but did 
not organize along the horizontal axis (unpublished data). Next 
we implemented kinesin-5, a well-characterized tetrameric, 
plus end-directed molecular motor essential for bipolar spindle 
organization in Xenopus (Fig. 1 C; G.Yang et al., 2007; Kapitein 
et al., 2008). Simulations with as few as five kinesin-5 motors 
per MT generated antiparallel arrays in which MTs underwent 
poleward sliding at rates near those measured for kinesin-5 mo- 
tility, while the motor itself was static in the spindle (Fig. 1 D; 
Videos 1 and 2). However, structures lacked any high density of 
aligned MT minus-ends that would be characteristic of spindle 
poles (Fig. 1 E). 

To generate poles, we modeled the effects of Xenopus 
kinesin-13s, kinesin superfamily protein 2a (Kif2a) and mitotic 
centromere-associated kinase (MCAK), which can depolymer- 
ize MTs from the minus end and are recruited to spindle poles 
by other pole-localized proteins in a dynein-dependent manner 
(Walczak et al., 1996; Desai et al, 1999; Gaetz and Kapoor, 
2004; lang et al., 2008). To localize the activity, we imple- 
mented the pole cohesion factor NuMA as an oligomeric MT 
cross-linker delivered to minus ends by dynein (Dionne et al., 
1999; Haren and Merdes, 2002). NuMA could bind MTs within 
0.5 um of a minus end (Fig. 2 A, yellow), or it could bind any- 
where along a microtubule, travel toward its minus end at the 
speed of dynein, and remain bound there (Fig. 2 A, cyan; Toba 
et al., 2006). Each NuMA protein recruited a fixed amount of 
kinesin- 13, which formed a gradient of depolymerization activ- 
ity around NuMA with maximal strength within a 0.5-um ra- 
dius, dropping to 10% of the maximal activity between 0.5 and 
1.5 um (Fig. 2 B, blue zone). As observed for kinesin-13s, the 
activity could bind to the MT lattice (Fig. 2 B, blue area of 
MTs), but could only depolymerize MTs from the minus end 
(Fig. 2 B, stars; Helenius et al., 2006). In the model, all kinesin- 13 
bound some portion of a MT, generating a minus-end depolymer- 
ization rate that depended on local MT density. The result 
was slow minus-end depolymerization in the midzone, where 
activity was low and spread over many MTs, and fast depoly- 
merization at the periphery, where activity was concentrated 
on a small number of MTs (Fig. 2 C). Simulations generated 
steady-state spindles with focused poles in which minus ends 
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Figure 1 . Self-assembly of a fluxing bipolar MT array with dynamic MTs and kinesin-5. (A] Simulations were initiated with MTs in an antiparallel array 
centered in the zone of high RanGTP concentration (yellow dashed line]. MT color indicates orientation (MTs with plus-ends pointed to the right in red, to 
the left in green) and region (minus-end region in white). (B) MTs experience repulsion at close range (0-50 nm) and attraction between 50 and 90 nm. The 
shaded area represents the interaction zone. (C) Homotetrameric kinesin-5 (orange) are simple plus end-directed motors, and slide apart antiparallel MTs. 

(D) Kymograph of selected MT minus-ends (green) and kinesin-5 motors (orange) shows that kinesin-5 remains stationary while bound MTs slide poleward. 

(E) Bipolar MT structure formed by MTs, cross-linking force, and kinesin-5. Single MTs extend far from the midzone, and poles are not defined. 



and NuMA converged sharply, hallmarks of spindle pole forma- 
tion (Fig. 2, C and D, blue and cyan lines, black arrows; Video 3; 
Merdes et al., 1996; Burbank et al., 2006). 

No model captures all features of the Xenopus meiotic spin- 
dle, and ours in particular lacks elements that determine spindle 
width. However, the model reproduced many key features of the 
organization of MTs along the longitudinal axis. The spindle 
length was well defined at 28 ± 0.86 urn (Table SI; Fig. SI A), 
and pole-to-pole MT organization appeared similar to Xenopus 
meiotic spindles (Heald et al., 1997; Z. Yang et al., 2007). Spindle 
length ranged from 24.4 to 32.0 urn when the number of MTs in 
the spindle was varied ± 50% (Fig. SI E), similar to experimental 
observations when the number of nonkinetochore MTs was re- 
duced in Xenopus extract (Houghtaling et al., 2009). Plots of MT 
density based on MT orientation (Fig. SIB) compared favorably 
to those interpolated from flux experiments (Burbank et al., 2006; 
G. Yang et al., 2007), as did fractional MT end density (Fig. SIC; 
Burbank et al., 2006). MT end density was also similar to experi- 
mental interpolations (Tirnauer et al., 2004; Burbank et al., 2006), 



with plus and minus ends enriched in the midzone and pole, re- 
spectively, but with a lower MT minus-end density at the poles in 
the simulation (Fig. SI D). This discrepancy suggests that nucle- 
ation near the spindle poles is higher in vivo than in the model. 

Kinesin-5 operated in the low load limit, near its maximum 
speed and similar to experimentally measured values (Table SI; 
Video 4). MT flux was spatially homogeneous (see discussion 
below) and occurred at experimentally observed rates (2.28 ± 
0.07 um/min [mean ± SEM]; Ohi et al., 2007; G. Yang et al., 

2007) . For simplicity, the simulated kinesin-5 did not have a 
preference for antiparallel MTs (van den Wildenberg et al., 2008), 
and was not transported by dynein (Kapoor and Mitchison, 2001; 
Uteng et al., 2008). The latter effect has been suggested to boost 
kinesin-5 bundling activity near the spindle poles (Uteng et al., 

2008) . This was not necessary in our model because the cross- 
linking force bundled MTs independently of the simulated 
kinesin-5-like motor. 

At spindle poles, the mean NuMA residency was 230 ± 26 s 
(SEM) with an average half life of 167 ± 30 s (SEM), similar to 
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Table I. Parameters of the simulation 



Parameter 

Total time 
Viscosity 

MTs 

Amount 
Steric repulsion 

Bundling attraction 

Plus-end assembly speeds 

Plus-end rates, low RanGTP 

Plus-end rates, high RanGTP 

Kinesin-5 motors 

Amount 
Motor domain 



Dual-motor complex 

NuMA sphere 

Amount 
Radius 

Turnover rate 

Binding to MT minus end 

Dynein driven motility 



Depolymerization of MT minus-ends 

Discretization parameters 

Time step 
MT section length 
MT initial length 
MT minimal length 



Value 

5,000 s 

0.1 pN s/um 2 

200 

Radius 0.025 urn 

Stiffness 500 pN/um every 0.25 pm 
Radius 0.045 urn 

Stiffness 50 pN/um every 0.25 pm 
Growth 0.18 pm/ s (1 0.8 pm/min) 

Shrinkage 0.33 pm/s (19.8 pm/min) 
Catastrophe 0. 1 s _1 

Rescue 0.004 s~' 
Catastrophe 0.024 s _1 

Rescue 0.01 8 s" 1 



6,000 

Max speed 0.04 pm/s (2.4 pm/min] 

Stall force 5 pN 

Binding range 0.05 pm 

Binding rate 2.5 s - ' 

Unbinding rate 0.05 s _1 
Stiffness 1 00 pN/pm 

Resting length 0.05 pm 

3,000 
0.025 pm 
0.1 s-' 

Max distance 0.5 pm 

Binding rate 1 s _1 
Max speed 0.8 pm/s (48 pm/min] 

Stall force 2.5 pN 

Binding range 0.05 pm 

Binding rate 0.5 s _1 

Unbinding rate 0.05 s _1 
—0.6 pm/min per sphere 

0.005 s 
0.5 pm 
0.1045 pm 
0.04 pm 



Note/ reference 

1 OOx water 

Kept constant by nucleation 
(Mitchison et al., 2004) 



(Carazo-Salas et al., 2001 ; 

Tirnauer et al., 2004) 
(Wilde etal., 2001; 

Brown et al., 2007) 
(Wilde etal., 2001) 



(Valentine et al., 2006; 
Korneevetal., 2007) 



Has no binding specificity 



Allows dissociation of aggregates 



(Toba etal., 2006) 



For bending elasticity 
Length when created 
Shorter MTs are deleted 



the half life of the mobile fraction of NuMA observed in cells 
(176 ± 46 s; Kisurina-Evgenieva et al., 2004). The simulation 
did not contain an immobile fraction of NuMA because NuMA 
de-oligomerization occurred at a constant rate. In agreement 
with in vitro observations of kinesin-13 activity on single MTs 
(Helenius et al., 2006), more than 99% of individual MT minus- 
end depolymerization rates were less than 4 um/min. As observed 
experimentally, flux rates were unchanged upon kinesin-13 in- 
hibition (2.22 ± 0.04 um/min [mean ± SEM]; Ohi et al., 2007). 

Inhibition of individual activities in the simulation re- 
produced experimental findings. For example, immobilizing a 
fraction of kinesin-5 in the simulation reproduced the inhibitory 
effects of the small molecule monastrol, slightly decreasing 
spindle length (Fig. S2 A; Burbank et al., 2007), and decreasing 
poleward flux rates (Fig. S2 C; Miyamoto et al., 2004) as long 
as the structure did not collapse. Keeping all other parameters 



intact but removing the ability of NuMA to oligomerize, or re- 
moving dynein-dependent transport of NuMA, disrupted spin- 
dle poles and lengthened structures, consistent with experimental 
results (Fig. S2, D-F; Merdes et al., 1996, 2000). 

Certain features, specifically those dependent on the MT- 
clustering ability of dynein, were not replicated by the model. 
In simulations, flux was spatially homogeneous (2.26 ± 0.06 
um/min [mean ± SEM]; unpublished data) and did not decrease 
near the spindle poles as seen in vivo (Burbank et al., 2007; 
Yang et al., 2008). At high monastrol concentrations, Xenopus 
extract spindles collapse in a dynein-dependent manner into a 
radially symmetric aster (Kapoor et al., 2000). In simulations as 
well, bipolar arrays collapse to form an astral structure (Fig. S2 A), 
although this pattern is not radially symmetric due to the resid- 
ual bundling interactions (Fig. S2 B). Furthermore, monastral 
spindle formation is mostly due to kinesin-13 activity. In the 
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Figure 2. Depolymerization recruited by dynein-transported NuMA generates spindle poles. (A) NuMA (yellow] is transported to MT minus-ends by dynein 
(cyan). Oligomerization of NuMA leads to its accumulation at the spindle poles. (B) NuAAA (yellow] recruits kinesin-1 3, whose binding (blue gradient) is 
spread among adjacent MTs. However, depolymerization (blue stars) only occurs at minus ends. Depolymerization is proportional to the amount of NuMA and 
inversely proportional to the local MT density (see Materials and methods). (C) Steady-state spindle structure with oligomerized NuMA (yellow) delivered to 
the pole by dynein (cyan). In the midzone, small amounts of kinesin-1 3 activity produce only slow depolymerization of minus ends. At the poles, accumulated 
NuMA recruits a large amount depolymerization activity distributed over few MTs, resulting in high depolymerization rates. (D) Minus-end density shows no 
polar accumulation (black) unless kinesin-1 3 depolymerization is added to the model (dark blue). With depolymerization, NuMA (light blue) accumulates at 
spindle poles (black arrows; average of 100, 95 simulations). (E) Schematic of metaphase steady state. MT density decreases from the midzone to the pole 
(green). Because depolymerization activity is spread over MTs, the depolymerization rate per MT increases with the distance from the midzone. 
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absence of kinesin-5-mediated flux, distal MT minus-ends are 
quickly depolymerized, causing NuMA oligomers to collapse 
into a single, central cluster, where dynein and depolymeriza- 
tion activities work together to align MT minus-ends. In addition, 
because the zone of primary MT nucleation and stabilization 
remains always centered in the simulations, the astral pattern 
more closely resembles the effects of kinesin-5 inhibition on 
chromatin bead spindles (Walczak et al., 1998). 

Mechanisms of spindle pole formation 

The establishment of a steady-state spindle requires polar MT 
disassembly to balance poleward MT flux (Fig. 2 E; Video 3). 
Such a balance is not obtained by fine-tuning the enzymatic 
rates of kinesin-5 and kinesin-13 (Goshima et al., 2005). Rather, 
Xenopus egg extract spindles exhibit robustness, forming over a 
range of flux and depolymerization rates, at the expense of 
changing spindle length (Shirasu-Hiza et al., 2004; Ohi et al., 
2007). A simple feedback mechanism confers a similar robust- 
ness in the simulation (Fig. S2 G, arrows): accumulated kinesin-13 
activity is spread over all MTs at the pole, so that the depoly- 
merization rate per MT depends on the number of MTs at the 
pole (Fig. 2 E). 

When the simulation begins, poles initially separate 
while depolymerization activity is spread over so many MTs 
that each MT slides faster than it depolymerizes (Fig. 2 E, top). 
As the poles extend farther, the number of MTs arriving at the 
pole decreases due to plus-end catastrophes. Eventually, the 
number of MTs arriving at the pole results in a depolymeriza- 
tion rate per MT that perfectly balances the flux rate (Fig. 2 E, 
middle). This feedback mechanism can also cause spindle 
poles to contract if MTs extend too far (Fig. 2 E, bottom), cre- 
ating a stable equilibrium for pole position (Fig. 2 E, middle; 
Video 5). Due to the stochasticity of MT dynamics and NuMA 
localization, spindle length fluctuates around the mean value 
(Fig. SI A). However, the kinesin-13 mechanism requires that 
MT density increase toward the midzone. In spindles com- 
posed of a constant number of MTs midzone-to-pole, such as 
in yeast, alternative feedback mechanisms must be used, such 
as MT length-dependent depolymerization by yeast kinesin-8 
(Varga et al., 2006). 

The MT depolymerization mechanism aligning minus 
ends is in principle independent from a motor-driven MT clus- 
tering (Burbank et al., 2006), yet both mechanisms seem to use 
the same motor in vivo: dynein. Theoretically, both mecha- 
nisms can work in concert, and neither alone fully explains 
spindle dynamics. Decreased rates of MT sliding have been 
measured for a "polar array" of MTs near the spindle pole in 
Xenopus spindles (Burbank et al., 2007; Yang et al., 2008), and 
arise in the slide-and-cluster model by antagonism between 
the dynein- and kinesin-5-like motors. The kinesin-13 mecha- 
nism does not replicate this phenomenon because dynein acts 
primarily to transport NuMA and is not abundant enough to 
oppose kinesin-5. On the other hand, the slide-and-cluster model 
incorporates only uncoordinated stochastic MT disappearance. 
In our model, flux is possible because minus-end depolymeriza- 
tion is linked to the arrival of MTs at the pole, where kinesin- 1 3 
is localized. 



MT amplification is essential for bipolarity 

How MT nucleation mechanisms contribute to spindle architec- 
ture is unclear. Nucleation within Xenopus spindles is not directly 
observable, as current reagents cannot distinguish between 
nucleation and plus-end rescue events, but -v-tubulin-mediated 
MT nucleation along existing MTs has been reported in Dro- 
sophila and human mitotic spindles (Goshima et al., 2008; Zhu 
et al., 2008). A similar MT-amplification mechanism in Xeno- 
pus would explain the exponential increase in MT mass during 
RanGTP-stimulated aster formation (Clausen and Ribbeck, 
2007). We therefore tested the sensitivity of the simulation to 
the percentage of RanGTP-mediated nucleation versus MT am- 
plification. We found the two sources of nucleation in the spin- 
dle to be highly redundant (Fig. SI F), likely reflecting the fact 
that the 20-um breadth of the RanGTP zone nearly matched that 
of MT amplification throughout the 28-um spindle. Restricting 
RanGTP-mediated nucleation from 20 to 10 um increased sen- 
sitivity to the nucleation mechanism, as structures required at 
least 75% MT amplification for bipolarity (Fig. 3 A, color indi- 
cates bipolarity). With greater than 25% RanGTP nucleation, 
MT structures were primarily astral, with minus ends close to the 
midzone, plus ends growing past minus ends (Fig. 3 B, bottom), 
and disrupted NuMA localization, illustrating the argument 
that plus-end growth disrupts bipolar MT organization if most 
nucleation occurs near the midzone. MT amplification also 
allowed large structures to achieve the "tiled array" composition 
proposed by G. Yang et al. (2007), in which short MTs exist 
throughout the spindle (Fig. 3 B, top). 

Spindles are robust to MT length 
distribution but sensitive to pole 
formation mechanism 

Estimates of the MT length distribution in Xenopus spindles 
vary significantly. Centrosomal assays have indicated an expo- 
nential length distribution dependent primarily on the plus-end 
dynamic instability parameters, whereas a statistical model 
from tubulin speckle measurements indicated a truncated nor- 
mal distribution within the spindle (Table SI). The simulation 
generates a truncated exponential distribution with a mean MT 
length of 6.51 ± 0.19 um (mean ± SE; Fig. SI G). 

When screening for sensitivity to different parameters, 
we noted that the mean MT length could change dramatically 
without disrupting spindle structure (Fig. 3, C and D), implying 
that spindle bipolarity is insensitive to shifts in the length distri- 
bution. However, spindles were highly sensitive to minus-end 
depolymerization, even though this mechanism accounted for 
only 3.2 ±0.1% (mean ± SE) of the total MT depolymerization 
in a simulation. Almost all MT depolymerization occurred at the 
plus end, as previously predicted (Ohi et al., 2007). Thus, the 
simulation demonstrates that a bipolar spindle could be com- 
posed of MTs with nearly exponentially distributed lengths. 

Length of simulated spindles scales with 
MT destabilization activity 

Next, we used the spindle simulation to examine which param- 
eters affected spindle length most. Varying the MT growth rate 
and nucleation mechanism did not have a strong effect (Fig. 3 D; 
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Figure 3. Spindle assembly is sensitive to mechanisms affecting internal MT organization but not to MT length distribution. (A) When the area of chromatin- 
mediated nucleation was decreased to span 10 pm, the bipolarity (green, yellow, and red) of MT structures was sensitive to the nucleation pathway. Color 
indicates MT "bipolarity index" defined in Materials and methods: astral (red), bipolar (yellow), strongly bipolar (green). MT structures remained bipolar 
only for percentages of chromatin-mediated nucleation up to 25%, while spindle length in the bipolar regime scaled slightly with nucleation pathway 
(92 simulations). (B) Position of the minus ends of right-oriented MTs at a single time point with a 1 0-um zone of chromatin-mediated MT nucleation (double 
arrow). In a simulation with 86% MT amplification (green), the bipolar MT structure contained short MTs throughout while long MTs were positioned with 
their minus ends at the pole (dashed lines). Four sample MTs are illustrated (arrows). In a simulation with only 25% MT amplification (red), the astral MT 
structure also contained small MTs throughout, but minus ends did not extend far from the midzone, and plus ends grew past minus ends. (C) Mean MT 
length and bipolarity index (color) for simulations in which parameters for the nucleation pathway, MT plus-end growth rate, MT plus-end catastrophe rate, 
and kinesin-13 rate were varied. Bold circles indicate the simulations represented in D (92, 99, 95, 99 simulations). (D) Total MT density along the pole- 
to-pole axis from representative simulations (bold circles in C), where color indicates bipolarity index. Spindles could form and fail in a variety of sizes, 
indicating a strong dependence for bipolarity on assembly mechanism rather than the exact organization of the MT structure. 



Fig. SI F), but spindle length varied dramatically with plus-end 
catastrophe frequency and minus-end depolymerization rate. 
Interestingly, these two parameters determine how far a MT 



could slide before disappearing. Bipolar spindles with lengths 
from 40 to 25 um formed within a relatively wide range of catas- 
trophe frequencies between 0.6 and 3.0 min -1 (Fig. 4, A and B), 
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structures obtained with different MT plus-end catastrophe frequencies (mean ± SD, 95 simulations). (C and D) Range of steady-state spindle structures 
obtained with different MT minus-end depolymerization rates (mean ± SD, 126 simulations). In A-D, the text or symbol color designates bipolarity index: 
astral (red), bipolar (yellow), strongly bipolar (green). 



which correspond to measured values (Wilde et al., 2001). 
Higher catastrophe frequencies led to aberrant structures with 
very short microtubules, whereas low frequencies were toler- 
ated due to the constant number of MTs in the simulation. Alter- 
ing the rate of MT minus-end depolymerization from 0.1 to 
1.2 um/min gave bipolar structures with an even wider range of 
lengths, from 45 to 17 um (Fig. 4, C and D), although simula- 
tions with very low activity did not form high MT density spin- 
dle poles. Interestingly, Xenopus spindle length has been 
documented to scale with titration of the kinesin-13s MCAK 
and Kif2a (Ohi et al., 2007) and the tubulin-sequestering pro- 
tein stathmin/Opl8 (Budde et al., 2001; Houghtaling et al., 
2009), which can act on both plus and minus ends to promote 
catastrophe or minus-end depolymerization. 

Conclusion 

The model discussed here offers a coherent picture of how MT 
dynamic instability, flux, and nucleation contribute to the self- 
organization of a steady-state bipolar MT structure. Simulations 
reconstituted a steady-state spindle composed of dynamic MTs 



using basic mechanisms: antiparallel MT sliding by a tetrameric 
motor, MT amplification, and minus-end cross-linking and de- 
polymerization. Because MT growth speeds exceed the speed 
of poleward motion, dispersed nucleation in the form of nucle- 
ation along existing MTs was necessary to organize dynamic 
MTs into a structure where MT plus-ends point inward. The 
simulation also illustrated that spindles can be composed of 
MTs with a truncated exponential length distribution, with short 
MTs present throughout the structure. Finally, whereas previous 
theoretical models proposed that spindle length is determined 
by a balance of forces between motors of opposite polarity, the 
length of simulated spindles in our model was determined by a 
balance between the speed of MT sliding associated with flux 
and the speed of MT depolymerization at the poles. 

The model makes several predictions that can be tested 
experimentally. For instance, simulations showed that plus- and 
minus-end dynamics independently influence spindle length, 
which may be verifiable once a specific inhibitor becomes avail- 
able, or more is discovered about the mechanisms of kinesin-13 
localization to the pole. Similarly, the contribution of MT 
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amplification to bipolar structure assembly can be explored 
when the Xenopus homologues of other known amplification 
proteins are characterized. Finally, the accuracy of our predicted 
spindle architecture and MT length distribution could be as- 
sessed using cryo-electron tomography and reconstruction of 
meiotic spindles. Although the model addresses pole-to-pole 
MT organization, construction of the spindle in the perpendicu- 
lar plane remains to be explored. Furthermore, the model could 
be extended to examine the contributions of additional ele- 
ments, such as chromatin and kinetochores, to the assembly of 
a functional spindle. 

Materials and methods 

Simulation 

We used Cytosim to solve the Langevin equation for the motion of flexible 
MTs and NUMA spheres at low Reynolds number (Nedelec and Foethke, 
2007]. Parameter values (Table I) were taken from experimental measure- 
ments when possible and otherwise broadly screened to identify regimes 
resulting in appropriate spindle morphology. Parameters varied include: 
strength of steric repulsion and attraction (not depicted), plus-end growth 
speed (Fig. 3, C and D), plus-end catastrophe (Fig. 3, C and D; Fig. 4, 
A and B) and rescue rates (not depicted), number of kinesin-5 motors (not 
depicted), speed/mobility of kinesin-5 motors (Fig. S2, A-C, G and H), pro- 
portions and area of nucleation mechanism (Fig. 3 A, Fig. SI F), amount 
of NuMA (not depicted), NuMA binding and transport efficiency (not de- 
picted), kinesin-13 depolymerization rates (Fig. 4, C and D), and MT and 
NuMA initialization conditions (not depicted). 

Simulations were run with a time step of 5 ms (similar results were 
obtained for smaller time steps), and iterations reached a total of 5,000 s 
unless otherwise noted. Measurements at steady state included only statis- 
tics from 2,500-5,000 s at 100-s intervals, except for kinesin-5 and dy- 
nein inhibition, where statistics were taken from 2,000-2,500 s. Spindle 
longitudinal axes were defined by the average MT orientation weighted by 
MT length. MT position and NuMA density were projected along this axis 
for line density measurements. Spindle length was calculated as the dis- 
tance between the points of maximal NuMA density on each side of the 
midzone (i.e., the spindle poles) with 1 urn resolution, and time-derived to 
yield the spindle elongation rate. Spindle width was calculated as the dis- 
tance spanned by MT ends at the midzone. MT density and MT end density 
were normalized by the number of MTs, and fractional MT end density was 
normalized by MT density. MT flux speed was calculated from the distance 
travelled along the longitudinal spindle axis by MTs during their lifetime 
and only included measurements for MTs with lifetimes >40 s. 

A bipolarity index in [—1,1] is defined from P(x) and A(x), the num- 
bers of MTs pointing to and away from the midzone, both functions of lon- 
gitudinal position x: 

B.l.= f dx\ \ 44 V4dx 

[is j } s P(x)+A(x) 

where 

S=[x\A(x)+P(x)>0\ 

is the segment covered by MTs. The bipolar index provides a quantitative 
measure of bipolarity for a MT structure: B.I. < 0 (red) indicates that 
MTs are predominantly oriented away from the midzone. On the contrary, 
B.I. > 0 for a spindle-like organization, in which MT minus ends are distant 
from the midzone while plus ends overlap in the midzone. We further differ- 
entiated between weak (0 < B.I. < 0.1 , yellow) and strong (0.1 < B.I., green) 
bipolar structures. 

MTs and cross-linking force 

MT plus-ends switched stochastically between growth and shrinkage, 
which occurred at constant speed for all MTs (Table I). The rates of switch- 
ing depended on the position of the plus end: either inside or outside the 
RanGTP zone (Table I). We obtained similar results with either a rectangu- 
lar or circular (Fig. SI E) RanGTP zone. MT minus-ends were stable unless 



depolymerized by the NuMA/kinesin-1 3 mechanism (see below). Any MT 
shorter than 0.04 pm was deleted from the simulation and replaced immedi- 
ately by a new growing MT of length ~0. 1 pm with a random orientation. 
The positions of new MTs were determined according to two nucleation 
pathways. In the first, MTs were nucleated at random positions within the 
RanGTP zone. In the second, MTs were nucleated at random locations on the 
side of existing microtubules. For each nucleation event, a random number 
was used to choose between the primary and secondary pathways, follow- 
ing their proportions (80% and 20%, respectively, except for Fig. 3 A). 

In the 2D simulations, MTs were not prevented from crossing, but a 
steric repulsive spring force with a resting length of 0.05 pm separated any 
MTs that came closer than this distance. The attractive cross-linking force 
also had a rest length of 0.05 pm but acted on any MTs closer than 0.09 pm. 
For repulsion and attraction, MTs interacted through points spaced at 0.25 pm 
along their length (Table I), with repulsion 1 0 times stronger than attraction. 
The attractive force represented approximately one cross-link of the strength 
of kinesin-5 per 0.5 pm of MT length. Forces produced by these interactions 
were projected perpendicular to the MT axis, so as to not generate any MT 
sliding forces, and thus did not impede kinesin-5-driven MT sliding. 

Motors 

Motor-MT attachment occurred at a constant rate, whereas the detachment 
rate increased exponentially with applied load: 

k 0 ff=koe 2///H "" 

Motors could bind and walk along MTs processively at a speed linearly 
dependent on load (Valentine et al., 2006; Gennerich et al., 2007): 

V=V max fl- / 1 



Kinesin-5 

Kinesin-5 was modeled as two separate motors connected by a linear spring 
of nonzero resting length, each walking toward MT plus-ends. Binding oc- 
curred independently of the orientation of MTs, and binding of the two motors 
to the same MT was forbidden. Unbinding was immediate upon reaching a 
MT end. The concentration of unbound kinesin-5 motors was considered to be 
uniform, and thus unbound diffusing motors were not simulated. 

Dynein and NuMA 

A single dynein was bound to each free NuMA and could also bind anywhere 
along MTs, traveling toward the minus end (Toba et al., 2006) and stalling at 
the final 0.5 pm of the MT. Because NuMA bound either a MT or a single dy- 
nein, only NuMA oligomers could cross-link MTs. Adjacent NuMA proteins 
could oligomerize via a Hookean force to form clusters with a 50-nm bond 
length. Turnover was implemented by displacing a fraction of the oligomerized 
NuMA to a random position in the simulation space at each time step. 

Minus-end depolymerization 

The depolymerization gradient was calculated at every time point by using 
a square grid size of 1 pm. The resolution of the grid was varied without 
affecting the simulation results. We first computed for each grid box a total 
amount of depolymerization activity proportional to the amount of NuMA 
within the region. Each NuMA recruited 0.01 pm/s depolymerization ac- 
tivity within 0.5 pm, and 0.001 pm/s activity within 0.5-1 .5 pm. Then, the 
depolymerization activity present in a grid box was shared equally among 
all MT segments present in the grid box. When a MT segment corre- 
sponded to the minus end of a MT, it depolymerized at the calculated rate. 
This procedure is inspired by the experimental observation that kinesin-13s 
can bind and diffuse along the MT lattice (Helenius et al., 2006). Activity 
from multiple NuMA proteins was additive, meaning that increased NuMA 
localization resulted in increased depolymerization rates. There is no 
mechanism to protect minus ends from this depolymerization activity. 

Initialization 

To minimize computation costs, simulations were initialized with horizontal 
microtubules (Fig. 1 A) and with 25% of MT plus-ends in the shrinking 
state. Simulations with random initial configurations eventually reached the 
same steady states (not depicted). NuMA could be initialized randomly 
throughout the simulation space or at the minus ends of microtubules, 
reaching similar steady states. 
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Online supplemental material 

Fig. SI presents simulated elements that are directly comparable with 
figures published elsewhere (see text for references), which contain mea- 
surements of unperturbed meiotic spindles. Fig. S2 illustrates standard 
perturbations such as inhibitions of essential molecular components of 
the spindle. Five videos display different aspects of the simulated spin- 
dles. A table compares spindle length, microtubule distribution, length, 
and flux speed, in the simulation and in Xenopus laevis and tropicalis. 
Online supplemental material is available at http://www.jcb.org/cgi/ 
content/full/jcb.201006076/DCl . 
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